arXiv:1507.08308v2 [physics.flu-dyn] 13 Dec 2015 


Mathematical model of influence of friction on 

the vortex motion 


Olga S. Rozanova a ’* Jui-Ling Yu b Chin-Kun Hu c 

a Department of Mechanics and Mathematrics, Moscow State University, Moscow 

119992 Russia 

b Department of Financial and Computational Mathematics, Providence 
University, Taichung, 49301, Taiwan 

and National Center of Theoretical Sciences at Taipei, Physics Division, National 
Taiwan University, Taipei 10611, Taiwan 

c Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan 
and National Center for Theoretical Sciences, National Tsing Hua University, 

Hsinchu 30013, Taiwan 


Abstract 

We study the influence of linear friction on the vortex motion in a non-viscous 
stratified compressible rotating media. Our method can be applied to describe the 
complex behavior of a tropical cyclone approaching land. In particular, we show 
that several features of the vortex in the atmosphere such as a significant track 
deflection, sudden decay and intensification, can be explained already by means of 
the simplest two dimensional barotropic model, which is a result of averaging over 
the height in the primitive equations of air motion in the atmosphere. Our theoretical 
considerations are in a good compliance with the experimental data. In contrast to 
other models, where first the additional physically reasonable simplifications are 
made, we deal with special solutions of the full system. Our method is able to 
explain the phenomenon of the cyclone attracting to the land and interaction of the 
cyclone with an island. 
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There exists a lot of theoretical, numerical, and experimental studies about the 
influence of land friction to the dynamics of tropical cyclones (e.g. [1] and ref- 
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erences therein). However, they contain sometimes contradictory results. For 
example, some theoretical and numerical studies of the sensitivity of tropic 
cyclones to the friction in axisymmetric models indicated that the intensity 
decreases markedly as the drag coefficient increases [4], [5]. At the same time, 
in [6], [7], [8], the authors present a series of three-dimensional convection- 
permitting numerical experiments in which the intensification rate and in¬ 
tensity of the vortex increase with the surface drag coefficient up to a cer¬ 
tain threshold value and then decrease. Further, the numerics made by MM5 
(Fifth-Generation Penn State/NCAR Mesoscale Model) have shown that the 
phenomenon of attraction of the cyclone to the land does exist [9]. Moreover, 
many experimental works indicate that the circular air motion in cyclone 
transforms into a convergent air flow during the landfall [11]. 

Of course, the structure of a tropical cyclone is three-dimensional and the pro¬ 
cesses of the moisture and heat transfer play an important role in its forma¬ 
tion. Nevertheless, we are going to show that several very complicated three- 
dimensional phenomena can be qualitatively explained by a relatively simple 
mechanical two-dimensional model where the topography is very crudely pa¬ 
rameterized by allowing the surface friction coefficient to vary. 

We deal with a special class of solutions of the full system, which is charac¬ 
terized by a linear profile of horizontal velocity. Under this condition we show 
that a smooth vortex in a material volume in three-dimensional stratified com¬ 
pressible non-barotropic flow and the simplest vortex in the two-dimensional 
barotropic model are governed by the same nonlinear system of ODEs. 

We are going to show that the following phenomena related to the cyclone 
motion can be explained by a rather simple low parametric model analytically: 
(i) a drastic intensification of the cyclone and a deflection of its track during 
landfall; (ii) complex behavior of the cyclone during interacting with the land 
(attraction and rounding); (iii) formation of converging flow starting from 
circular motion inside the cyclone. 


1 Full model of gas dynamics adapted for atmosphere 


We consider the system of non-isentropic polytropic gas dynamics equations 
in a uniformly rotating reference frame for unknown functions p > 0, p > 0, 
U = (Ui,U- 2 ,Us), S (density, pressure, velocity and entropy), in the presence 
of a horizontal dry friction, namely 

p(d t U + (■ U , V)U + le 3 xU + pU H + ge 3 ) = -Vp, (1.1) 

d t p + div(pU) = 0, (1.2) 
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d t s + (U, VS) = 0 . (1.3) 

The functions depend on time t and on point x € M 3 , e 3 = (0,0,1) is the 
’’upward” unit vector, l is the Coriolis parameter, p is the friction coefficient, 
g is the acceleration due to gravity (points in —e 3 direction), Uh = (fd, U 2 , 0), 
p > 0 is the friction coefficient. The state equation is 

P = d 7 e 5 , (1.4) 

where 7 > 1 is the adiabatic exponent. This system is traditional, see e.g. [12], 

[13]- 

For onr convenience taking into account (1.4) we write the equation (1.3) in 
terms of pressure: 

d t p + (U , Vp ) + 7 pdwU = 0. (1.5) 


Let us consider classical solutions of (1.1), (1-2), (1.5). For p = 0, the system 
implies the conservation of mass A4 = f pdx, momentum V — f pU dx 

£l(t) Q(t) 


and energy 



dx, 


inside a material volume O(t), if we assume the hydrostatic balance 


= -gp- 


( 1 . 6 ) 


To prove these conservation laws we apply the formula for the derivative with 
respect to time of integral taken over a material volume [14], namely, 

J f(t, x)dx — I x) + div(, f(t,x)U)) dx. (1.7) 

£l(t) fi(t) 


Let us introduce the following functionals: 

G ( t , x 3) = ^ I p\Xi\ 2 dxidx 2 , 

n H (t) 

G Xl (t,x 3 ) = l > f px\dx 1 dx 2 , 
fMb 

G XlX2 (d *£ 3 ) = g 

where Ad = (xi,x 2 ), X 2 = (x 2 ,— X\), i = 1,2, fl^(f) is a section of f2(t) at 
a fixed x 3 . We note that G{t,x 3 ) > 0 and A {t,x 3 ) = G X1 G X2 — G 2 lX2 > 0 for 
nontrivial solutions to (1.1), (1-3), (1.5). 


Fi(t,x 3 )= J (U,X i )pdx 1 dx 2 , 

(t) 

G X2 (t,x 3 ) = ^ f px 2 2 dxidx 2 , 
n H (t) 

J px\x 2 dx\dx 2 , 
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Let us consider l = const > 0 and // = p,(x 3 ). Further, we assume U 3 = 
ty(t,x 3 ). The last assumption, (1.1), and (1.6) imply that U 3 (t,x 3 ) satishes 
the Hopf equation 

d t U 3 + U 3 d X3 U 3 = 0, x 3 > 0. 

It is well known that if the initial data U 3 (0,x 3 ) are not increasing, then the 
solution will lose the smoothness within a finite time. However, the natural 
conditions U 3 (t, 0) = U 3 (t,+ 00 ) = 0 imply that the initial datum U 3 (0,x 3 ) is 
not increasing everywhere. Therefore the only possibility to consider a smooth 
solution for all t > 0 is to set U 3 = 0. This is the asymptotic of the solution 
containing shock waves as t — * 00 on ( 0 ,+ 00 ) for any initial data f/ 3 ( 0 ,a; 3 ) 
with the property t/ 3 (0,0) = /7 3 (0,+oo) = 0. Thus, the assumption U 3 — 0 
satishes this limit case. 

Lemma 1.1 For the classical solutions of (1.1), (1-2), (1.5) with U 3 = 0 the 
following relations hold: 


d t G — Fi, d t F 2 = lF l ~pF 2 , 


dfF\ = 2(7 — 1 )E p + 2 E k — IF 2 — ///‘V§ 
dfE = —2 pE k , 

where E k (t,x 3 ) = 1 / p\U h \ 2 dxidx 2 , E p (t,x 3 ) = -d— j pdxidx 2 , E = 
n H (t) 7 n H (t) 

E k + E p . 


Proof. To prove the identities it is enough to apply formula (1.7) with respect 
to the variables x 2 and x 2 . For example, taking into account (1.2), we get 


d t G = * I d t p\Xi \ 2 dx\dx 2 = 
n H (t) 

J ^ dw(pU) \Xi \ 2 + ^ div (p Uh |Ad| 2 ))^) dxidx 2 = 
vlhF) 

= j (X] , U) pdx\dx 2 = F\. 

The proof of other identities are analogous. □ 

Let us assume a special structure of velocity inside fl(t). Namely, we set 


Uh — QhXi, Qh 


( a H {t,x 3 ) b H (t,x 3 ) 
y c H (t, x 3 ) d H (t,x 3 ) 


U 3 = 0. 


( 1 . 8 ) 


It is well known from experimental data that the profile of velocity near the 
center of atmospherical vortex like tropical cyclone is approximately linear 
[15]. Thus, we can consider the core of the cyclone as Q(t). 
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Lemma 1.2 For the velocity ( 1 . 8 ) we have 


dtG .ri 2 q>hG Xi + 2 huG x ^x2t d)tG X2 2duG X2 T 2 chG X iX2: 

dtG ' xlX2 = («f/ + ( 1 h)G X\X 2 + huG X2 + chG Xi 5 

= —(7 — 1 )(oj? + dH)E p , d t A = 2 (a H + < 7 f/)A. 
The potential energy E p is connected with A as 

E p (t,x 3) = £ p ( 0 , i 3 )A (7 “ 1)/2 ( 0 , i 3 ) A^ 7+1)/2 (i, i 3 ). 


Proof. The proof is a direct computation with taking into account formula 
( 1 . 7 ). The variable 2:3 plays a role of parameter. □ 

Let us introduce new functions 

Gi = G X1 A~ (7+1)/2 , Ct2 = Gx2 A-l^ 1 )/ 2 , G 3 = G X1X2 A~G+P/\ 

Lemmas 1.1 and 1.2 imply that for the elements of the matrix Q and G 1, G 2 , G 3 
can be obtained the following closed system of equations: 


d t G\ 

= ((! ~l)a H ~ 

(1 + 7 )dn)G] 

+ 2 bjjG[ 

d t G 2 

= (( 1 - 7 )d H ~ 

(1 + 7 )a H )Gr 

> + 2 c h G- 

d t G 3 

— chG\ + bnG 2 

- 7 ( a H 

+ d.H 

)G?„ 

dtUrt 

= —a 2 H — b H cn + Ich ~ 

h a H ' 

- 1CG 2 , 

dtf>H 

= bn {p>H + d H 

) + Id.H 

— pb H + ICG 

d t c H 

= —ch(o,h + d.H 

) — lO-H 

— p.cn + ICG 3, 

dtdn 

= — d 2 H — b H cu - 

- Ibn — 

/idn ■ 

- ICG 1 , 

with /C = 

F^E p AG-V/\ 

=0- 




Thus, ( 1 . 9 ) describes the behavior of full 3 D system of dynamics of strati¬ 
fied atmosphere near the center of atmospherical vortex.For every fixed level 
£3 = 2:3 we have its own motion of air and the dynamics of the whole ma¬ 
terial volume is determined by a(0,a;3), 6(0, x 3 ), €(0,2:3), <^(0,2:3), (0,2:3), 

G 2 ( 0 ,2:3), G 3 ( 0 ,2:3). 


Below we are going to show that the same system of ODE describes a behav¬ 
ior of the simplest possible solution of barotropic two-dimensional model of 
atmosphere. 
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2 Bidimensional models of the atmosphere 


Since the horizontal scale in the atmospherical motion is much larger than the 
vertical scale, there exist many approaches to simplifying the model [16], [17]. 
Moreover, there is a possibility of averaging over the height to hide vertical pro¬ 
cesses and reduce the primitive system of equations to two space coordinates 
(see [18] for barotropic case and [19] for general case). Let us recall shortly 
the procedure of averaging. Let p, U = (Pi, U 2 , U 3 ), p, denote in the three- 
dimensional density, velocity and pressure. Namely, all these functions depend 
on (t,xi,x 2 ,x 3 ), x 3 G K+. Let us introduce 0 and / to represent an average 
of 0 and / over the height, respectively. The averaged values are introduced 

/s r°° — 1 /•oo 

as follows: 0 := / (j)dx 3 , f — pf dx 3 , where 0 and / are arbitrary 

Jo p Jo 

functions, and denote g(t,xi,x 2 ) = p, P(t,x\,x 2 ) =0, U(f, £ 1 , 0 : 2 ) = (Pi, P 2 ). 

Moreover, the usual adiabatic exponent, 7, is related to the “two-dimensional” 

27 -I 

adiabatic exponent 7 h as follows: 7 h = - < 7. 

7 

The impenetrability conditions are included in the model. We require the 
vanishing of derivatives of the velocity on the Earth surface and a sufficiently 
rapid decay for all thermodynamic quantities as the vertical coordinate x 3 
approaches to infinity. In other words, the impenetrability conditions assure 
the boundedness of the mass, energy, and momentum in the air column. They 
also provide the necessary conditions for the convergence of integrals. 

If / and /i are constant, the resulting two-dimensional system consists of three 
equations for density g(t,x), velocity U(f,x) and pressure P(t, x), x e M 2 : 


g(d t U + (U • V)U + £U) + VP = 0, 

( 2 . 1 ) 

d t g + div(pU) = 0, 

( 2 . 2 ) 

d t P + (U ■ VP) + j H P divU = 0. 

(2.3) 


Here C = IL + p,I, L 



, / is the identity matrix, 7 # G ( 1 , 2 ) . 


We used this model with p, — 0 in our previous papers [20], [21]. 
We can apply to (2.1)-(2.3) all considerations of Sec.l. 
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However, for us it will be convenient to restrict ourselves to the barotropic 
case, where P = Cp 7 , C = const. Thus, the system under consideration can be 
reduced to two equations (2.1), (2.2). 


7 — 1 

We introduce a new variable n = P -> and get the system 

d t U+ (U- V)U + (IL +/x/)u + c 0 V7r = 0, 
d t ir + (V7T • U) + (7 — 1) 7r div U = 0, 


(2.4) 


1 

with c 0 = y . 


2.1 A class of exact solutions 


We consider a simple class of exact solutions which correspond to the first 
terms of expansion of n at a critical point. Namely, we look for the solution 
in the form 


U(t,x) = Qx, 


Q 


( a (t) 

\c(t) d{t)) 


(2.5) 


7T(t, x) = A{t)x\ + B(t)x 1X2 + C{t)x 2, 

to get a closed ODE system for the components of the matrices Q and R = 

( A(t) |B(t)Y 
\kB{t) C{t),J 


R + RQ + Q t R + (7 — l)tr QR = 0, 
Q + Q 2 + ILQ + pQ + 2cqR = 0. 


( 2 . 6 ) 


The system of matrix equations consists of 7 nonlinear ODEs and has a very 
complicated behavior. First of all, we discuss the behavior of the system in 
the frictionless case. In fact, an analogous class of solutions was considered in 
[23] in another context. 

It can be readily checked that system (2.6), coincides formally with (1.9) at a 
fixed x-,i , where Q = Qh, 2co = /C, A = G 2, B = — 2G3, C = G\. Nevertheless, 
the nature of these systems is different and the components of solution of (2.6) 
are not a result of integration of corresponding components of solution of (1.9) 
with respect to x 3 . In particular, the constants 7 and 7 h are different for these 
systems. 
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2.2 A friction-free vortex (/i — 0), stability issue 


2.2.1 Axisymmetric case [20] 


It is easy to see that (2.6) has a closed submanifold of solutions having addi¬ 
tional properties a = d, c = —6, A — C, B — 0. These solutions corresponds 
to the axisymmetric motion. Note that it is the most interesting case related 
to the vortex in atmosphere. Here we get a system of 3 ODEs: 

A + 2^aA = 0, 

a + a 2 — b 2 + lb + 2c 0 A = 0. (2.7) 

b + 2 ab — la = 0, 


The functions a, b, A > 0 correspond to one half of divergence, one half of 
vorticity and the fall of pressure in the center of vortex respectively. The only 
nontrivial equilibrium point that relates to a vortex motion is 


a = 0, b = —c = b*, A = A* 


b*(b* — /) 
2cq 


The center of vortex corresponds to a domain of low pressure only if A* > 0 
(the motion is cyclonic). This implies b* < 0 or b* > l. Further, there exists 
one first integral 

b= l - + kA t, (2.8) 

where k is a constant [20]. Thus, (2.7) can be reduced to the following system: 


A = — 2'yaA, 


l 2 2 

a = — a 2 — — + k 2 A~i — 2c 0 A. 


On the phase plane {(H,a), A > 0}, there always exists a unique equilibrium 
(A*, a*) = (H 0 ,0), stable in the Lyapunov sense (a center), where A 0 is a 
positive root of equation 

l 2 9 2 

— + 2c 0 A = k 2 A-y. 


2.2.2 General case 

As follows from [24], the axisymmetric form of 2D vortex is stable with re¬ 
spect to asymmetric perturbations for the solution to the incompressible Euler 
equations. Indeed, the incompressibility condition implies aft ) + dft) = 0 and 
this reduces the full system (2.6) to (2.7). As we have shown in Sec.2.2.1, the 
equilibrium in this case is stable for any b* and l. 

Nevertheless, in the compressible case this property does not hold for arbitrary 
values of parameters. 



Theorem 2.1 [22] If 


. * l-V2 r ,* 1 + V2, , 

6 < - l or b > - l > 1, 

2 2 


£/ien £/ie equilibrium of system (2.6) is unstable. 

Proof. As one can readily check, the point 

a = cl = 0, b = -c = b*, A = C = A* = - ^ , 5 = 0 

2co 

is the only equilibrium of the full system (2.6). It is the same point of equilib¬ 
rium as in the axisymmetric case (2.7). Nevertheless, in the symmetric case 
this equilibrium is always stable in the Lyapunov sense, whereas in the general 
case the situation is different. Indeed, the eigenvalues of matrix corresponding 
to the linearization at the equilibrium point are the following: 


Ai = 0, A 2 , 3 = ±V~(2(2 - 7 Mb* ~l) + l 2 ), 


A4,5,6,7 — 

Since (2 — y)6*(6* — /) + l 2 > 0 for 7 e (1,2), then 9?(A 2 ,3) = 0. Eigenvalues 
A^ i = 4, 5, 6, 7 have zero real part if and only if b* satisfies the following 

inequalities simultaneously: l(b* + |) > 0, j + b*l — ( b*) 2 > 0, l 2 (b* + £) > 

(b* + (y + b*l — (&*) 2 ) , that is b* G 1 b/ 2 /, 1+ 2 ^ 2 l . For others values of 
b* the eigenvalues A 4) 5 j6) 7 = ±a ± i/3, a ^ 0, /3 ^ 0, therefore there exist 
an eigenvalue with a positive real part. Thus, the Lyapunov theorem implies 
instability of the equilibrium for b* < l and b* > l > l. □ 

Remark 2.1 We notice that the full system (2.6) has the first integral 
b — c — l — const V = A AC — B 2 . 

This reduces (2.6) to the system of 6 equations. If b* E E, E = o) U 

/, 1+ /' 2 l'j, then the matrix, corresponding to the system, linearized at the equi¬ 
librium, has 3 pairs of pure imaginary complex conjugate roots i — 2,..., 7 
(for the range of parameters under consideration the roots are simple). A study 
of stability in this case is extremely complicated and we will not dwell here (see 
[ 22 ]). 



9 



2.3 Influence of the friction on an axisymmetric vortex 


The system of equations describing a vortex with a rotational symmetry is the 
following: 

A + 2 r yaA = 0, 

a T a 2 — b 2 T lb T 2c 0 A = — pci., (2.9) 

b + 2 ab — la = —pb, 

The solution to the equation has a complicated behavior, nevertheless it is 
possible to study it analytically to a certain extent. 

Theorem 2.2 System (2.9) has two equilibriums (a\,b\,A\) = (0,0,0) and 
{a* 2 ,b* 2 , A* 2 ) = (—p,l, 0), both are unstable. 

Proof. Indeed, the matrix of the system linearized at the point (A 0 , do, b 0 ) is 

^-2"u„ -2 jAq 0 

Q(A 0 ,a 0 ,b 0 ) = — 2 cq —2ao~p 2b$ — l 

^ 0 — 2bf) T l —2a 0 — p 

The eigenvalues of Q( 0, — p,l ) solve the equation 

R(k) = (2j p — k)((p — k ) 2 + l 2 ) = 0. 

The polynomial has a positive root. This means instability of equilibrium 
(— p, /, 0). The eigenvalues of Q( 0, 0,0) are (0, — p±il ), therefore the linearized 
theory does not give an answer to the question about the stability or instability 
of zero equilibrium. However in the critical case we can use the theory of [25], 
Sec.4. Namely, we consider expansions into series a(A) = a±A + 0(A 2 ) and 
b(A) = biA + 0(A 2 ) as A — » 0. Then we substitute the expansions into (2.9) 
and get a 1 = h = -fpfp. Therefore A = A 2 + 0(A 2 ). This 

implies instability of zero equilibrium. □ 

Remark 2.2 Theorem 2.2 implies that the the zero equilibrium of the full 
system (2.6) is also unstable. 

Theorem 2.3 The solutions to system (2.9) has no finite time blow up points 
at t > 0 and the following inequalities hold: 

A < — C °—A + k 0 A^e~ 2 ^\ (2.10) 

7-1 

A < Rfe^, (2.11) 

where A = a2 f b 2 , S 0 , K 0 and k 0 are positive constants depending only on initial 
data. 
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Proof. First of all we note that the first equation and two latter equations of 
(2.9) imply 

A~ lh - 2aA- 1 h = 0 (2.12) 

and 

A + 2aA + 2c 0 aA + 2/xA = 0, (2.13) 

respectively. Equations (2.13) and (2.12) result 

= -2//AA -2 < 0. (2.14) 

From (2.14) we obtain (2.10). 


Let us prove inequality (2.11). First of all we note that inequality (2.10) implies 
that there exists a constant A, depending on initial data such that for A > A 
we have 

A < k x A (2.15) 

with a positive constant k\. Let us introduce a new variable W = {b — 
k)A~^e tlt . It is easy to check that 


dW 

dt 


< 0 . 

2 


(2.16) 


As follows from (2.16), Wft) < W( 0). The second equation of (2.9) takes the 
form 

a' = -a 2 — /ra + w 2 A 2 ^e~ 2 ' rflt - 2c 0 A - l ~. (2.17) 

First, we consider the cases In = 0, where Wft) = W( 0), and 6(0) < | (or 
W( 0) < 0), for In A 0- Then (2.17) and (2.15) imply that for sufficiently large 
A we have 

a' > fE 2 (0)A 2/7 e- 2wt - k 2 A, k 2 > 0. (2.18) 

Thus, if there exits an interval of t such that the inequality 

A^T > k 3 e 2 ^\ k 3 = k 2 /W 2 { 0), (2.19) 


then for these t the function aft) increases and, as follows from the first equa¬ 
tion of (2.9), Aft) decreases. Thus Aft) increases if and only if inequality 
(2.11), opposite to (2.19), holds. 

The last case is In A 0, 6(0) > | or W(0) > 0. Due to (2.16) there are two 
possibilities: Wft*) < 0 for some 6* > 0 or Wft) > Wo = const > 0 for all 
t > 0. The first case can be reduced to the case W(0) < 0 if we take 6* as 
the initial moment of time. We note that W ft) cannot be identically zero, this 
contradicts to the last equation of (2.9). If we assume W 0 > 0, the second 
possibility implies inequality (2.18) with Wf instead of W 2 ( 0). Thus, estimate 
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(2.11) follows from the same reasoning. Let us show that Wq does not vanish 
and our assumption is correct. Indeed, from (2.15) we have 


W < k^A^-y e^, k/± = const > 0, 


( 2 . 20 ) 


for sufficiently large A. Further, from (2.16) and (2.20) we obtain 

dXV 2 7j jt , 

> —k§W 2 -^e 2 —<, fc.5 = const > 0. 


dt 


( 2 . 21 ) 


We divide variables in (2.21) and integrate. After obvious estimates we get 

W{ 0) 


W > 




:= Wo > 0 . 


Thus, (2.11) is proved. □ 


Remark 2.3 For p = 0 inequalities (2.11) and (2.10) imply that the solution 
to system (2.7) is bounded for all t > 0 by a constant depending on the initial 
data. 


Remark 2.4 As follows from (2.16), the value ofW(t) is constant for l p = 0. 
From the conservation ofW for p = 0 we get the integral (2.8). 

Remark 2.5 Although systems (2.6) and (1.9) at a fixed x 3 are formally 
equivalent, we can obtain from (1.9) more information. Indeed, (2.6) is con¬ 
sidered in the whole space M 2 and it is not based on conservation laws. In 
contrast, (1.9) is considered in a moving volume, where the balance of energy 
E'(t ) = —2 pE k ft) holds. Let us recall that the component Aft) of solution of 
(2.6) corresponds to G\ in the solution of (1.9). Since we deal with the ax- 
isymmetric case, G\ = G 2 = {G/ 2) -7 . Let us note that if the velocity field has 
the form (2.5), a = d, c = —6, then E k = (a 2 + b 2 )G > 0, E p = diG 1 ^ 1 > 0, 
E k + E p < 62 , where and 62 are positive constants, depending on initial data. 
This implies G > (^i/^)^ 1 := £3 > 0. Thus, G\ = (G/ 2) -7 < (5 3 /2) -7 . In 
terms of system (2.6) this means that A is bounded. 

Figs.l shows the behavior of functions Aft), related to the intensity for dif¬ 
ferent periods of time (0-2 days, 2-2.3 days). The parameters are typical 
for geophysical vortex near its center. Namely, c 0 = 0.1, 7 = |, / = 7.3 • 
10 -5 s -1 , 6(0) = —5 • 10 -6 s -1 , A(0) = 1.95 • 10~ 7 (solution to equation 6 2 (0) — 
/6(0) — 2c 0 A(0) = 0). The parameter p is 1 • 10 -4 s -1 . The initial data cor¬ 
respond to a steady vortex for p = 0. One can see that in the presence of 
the land friction, the vortex first intensifies, after that a very fast oscillations 
begin. These oscillations can be interpreted as a destruction of the vortex. 
Moreover, the increasing amplitude of oscillations contradicts to boundedness 
of A in a moving volume (Remark 2.5). Nevertheless, this contradicts only 


12 



the smoothness of solution in the moving volume and implies the formation 
of frontal zone inside it. 

Further, Figs.2 show the field of wind for a steady vortex in the non-frictional 
case and the respective field influenced by the constant surface friction for 
the period of intensification of vortex, as in the left Fig.l. We can see a for¬ 
mation of a convergent stream. Fig.2, right, is taken from [11] and presents 
the experimental evidence of the fact that the streamlines of the wind in a 
tropical cyclone form the focal point near the landfall. Thus, in the frame of 
the model it is possible to explain the formation of a convergent stream from 
an axisymmetric steady state. 

Remark 2.6 The convergent inflow into cyclone is a well known feature of 
rotating boundary layers in general and can be explained be increasing of the 
Ekman pumping. Nevertheless this phenomenon can be explained within a two- 
dimensional non-viscous model. 




Figure 1. The intensification (left: 0-2 days) and destruction (right: 2-2.3 days) 
of the vortex, A(t) (intensity), the time is in seconds. 
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Figure 2. A convergent stream influenced by friction. Left: field of wind (2.5), ax¬ 
isymmetric case, built for a(t) and b(t) (see(2.7)), the period of intensification. Right: 
streamline analysis for the hurricane Frederic, experimental data [11] 
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3 Local and leading fields separation 


Let us change the coordinate system of (2.4) in such a way that the origin of 
the new system x = (xi,x 2 ) is located at a point X(t) = (Xi(t),X 2 (t)). It 
is associated with the center of vortex (here and below we use the lowercase 
letters for x to denote the local coordinate system). Let U = u + V, where 
v (t) = (yi(t),y 2 (t)) = (*!(*),x 2 (t)). Thus, 

d t u + (u • V)u+ 

"V + (/ L + /i/)(u + V) + coV7r = 0, 

d t TT + (V7T • u) + (7 — 1) 7T div U = 0. 

Given a vector V, the trajectory of vortex can be found by integrating the 
system 

x\{t) = V 1 (t), X 2 (t) = V 2 (t). (3.1) 

We assume that the pressure field can be separated into two parts 

7T = TT 0 (t,Xi,X 2 ) + 7Ti(t,X 1 ,X 2 ). 

The first field (we will call it local) is associated with the vortex, the second 
field can be considered as a leading one. In fact, the local field can be con¬ 
sidered as a perturbation of the leading field due to the vortex. We impose a 
requirement 

V 7T 0 | x =0 = 0. 

Formally we can write 

d t u + (u • V)u + £u + c 0 V 7 r 0 + 

V + C V + C 0 V7Tij = 0, 

r 1 (3-2) 

d t 7 T 0 + (V 7 T 0 • u) + (7 - 1) 7 To div U + 
d t TTi + (V 7 Ti • u) + (7 — 1 ) 7 Ti clivu = 0, 

Let us denote 

q = c 0 V 7 T!(f,x) - V 7 Ti(f,x) . (3.3) 

L x=oJ 

If we solve separately the system for the local field 

d t u + (u ■ V)u + £. u + c 0 V 7 T 0 + q = 0, 
d t n 0 + (V7T 0 • u) + (7 - 1) 7r 0 div u = 0, 

we get a linear equation for 717 : 

d t Ti 1 + (V 7 Ti • u) + 7 Tidivu = 0, (3.4) 
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which can be solved for any initial condition 7Ti(0, x). 
Further, from (3.2) we obtain 


V (t) + £V (t) + c 0 V 7 r 1 (i, x) 


x=0 


0 , 


then (3.5) and (3.3) result 


(3.5) 


X(t) + £X(t) + C 0 V7Ti(t, x) 


x=0 


0 . 


If q — 0, then we obtain a complete separation of two processes and the local 
held does not depend of the leading held. It is easy to see that q — 0 if and 
only if Hi is linear with respect to the space variables. If q ^ 0, however it is in 
some sense small, we can talk about an approximate separation of processes; 
\q\ plays a role of measure of separability of the local and leading processes. 


4 Influence of friction on the trajectory of vortex 


If the couple (u, 7To) is found, then the position of the vortex can be deter¬ 
mined from linear equations (3.4) and (3.3). Let (u, 7 r 0 ) be the exact solution 
considered in Sec. 2 . 1 . As follows from (3.4), if 

7Ti(0, xi,x 2 ) = M(0)xi + N( 0)x 2 + K( 0), 

then the initial held Hi is linear with respect to the space variables, that is 

!Ti(t, Xi,x 2 ) = M(t)x 1 + N{t)x 2 + K(t). 

Thus, we obtain the zero discrepancy term q (i.e. the complete separation of 
local and leading Acids). The coefficients M(t), N(t ) and K(t ) can be found 
from the ODE system 

M + (2y — 1 )aM — bN = 0, 

N + (2y — l)aN + bM — 0, (4.1) 

K + 2(7 — l)tr QK = 0. 

To obtain the trajectory of the center of vortex (Ad(t), X 2 (t)), we have to solve 
the system (3.5), which can be reduced to 

V\ — IV 2 + fiV j + CqM = 0, 

V 2 + IV j + /iV, + CqN = 0. 

Then the trajectory can be found from (3.1). The coefficients M(t ) and N{t) 
can be considered as a measure of intensity of the leading field. It is natural 
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to assume that they are so small that the local held can be discerned in the 
leading held (see the numerical examples from [21]). The function K(t) does 
not influence the trajectory. 


4-0.1 Constant coefficient of surface friction 

As we have shown, the friction basically causes the intensification of vortex. 
At the same time, the trajectory of vortex changes according to the lead¬ 
ing held and initial velocity. It can shrink or amplify, formation of loops is 
a typical behavior. Fig.3 shows the position of vortex, computed for the pa¬ 
rameters corresponding to a tropical cyclone within two days for /jl = 0 and 
H — 2 • 1CT 5 s -1 respectively. The Coriolis parameter l = 7.3 x 10 -5 s -1 , that 
corresponds to the latitude 30° approximately, c 0 = 0.1 (appropriate dimen¬ 
sion), 7=| (recall that in the procedure of averaging over the height, the 
value of heat ratio for air changes). Initial data are the same for the both 
cases, namely, Vi(0) = 14(0) = 1 m/s, M( 0) = N(0) = 10 -3 (appropriate 
dimension), b* = —10 _6 s _1 , initial condition corresponds to equilibrium for 

H — 0 . 



Figure 3. Influence of friction on the trajectories of ’’cyclone” within 2 days: 1 - 
H = 0, 2 - /t = 2 • 10~ 5 m/s, initial data are identical. 

4-0.2 Influence of a land: interaction of vortex with the ’’island” 

The most intriguing phenomenon is the interaction of the cyclone with the 
land. There are a lot of experimental and numerical evidences showing that 
the cyclone ”feels” the land; it can be attracted to the shore, but some¬ 
times it ’’avoids” the shore, on the contrary [10], [2], [3]. Now we are go¬ 
ing to show that the approximation of trajectory made by the ODE system 
(2.9),(4.1),(4.2), (3.1) can reproduce this complicated behavior. The coefficient 
of surface friction /i in this experiments is a function of space variables X \, X 2 . 
It changes from zero (sea surface) to some constant value /i 0 (’’island”). Let 
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li(X u X 2 ) — u(Xi,X2)fJ,o, 


u{X u X 2 ) = ^arctan _ arctan ^ 

where a is a very small constant. The square [—ah, ah] x [—£ 2 , £ 2 ] corresponds 
to the island, the motion begins over the sea. 

Fig.4 shows that the attraction to the lands increases with the roughness of 
the land. Fig.5 shows the difference of divergence and intensity for different 
roughness of the island. Fig.6 shows that for some initial conditions, the cy¬ 
clone can avoid the island. 




Figure 4. Attraction of cyclone to the island depends on the coefficient of friction. 
Computations made for 5 days, Li(0) = 1^(0) = lrn/s, M(0) = N(0 ) = 10 -3 
(appropriate dimension) Left: 1 - po = 0, 2 - po = 10~ 4 , 3 - po = 2 (s _1 ). 




Figure 5. Intensity (left) and divergence (right) of the vortex. 1 - po = 10 5 , 2 - 
po = 2 • 10" 5 , 3 - po = 2.5 • 10~ 5 (s _1 ). 


5 Conclusion 


To study different properties of vortices in a compressible medium (e.g. at¬ 
mosphere), we analyze a special class of smooth motions characterized by a 
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Figure 6. Avoiding of the shore: comparison with the motion over sea, dimen¬ 
sions of the island are [—70,70] X [—200,200] (km). V] (0) = V^O) = lm/s, 
M( 0) = —10 -4 , IV (0) = 10” 3 (appropriate dimension), computations made for 
8 days. 1 - /xo = 0, 2 - /j,q = 10 -3 (s _1 ). 

linear profile of the horizontal velocity. It is well known that the velocity has 
this property near the center of vortex (e.g.[15]). We consider two models: 
the primitive 3D model of atmosphere and 2D model obtained from it by the 
standard averaging procedure. 

In the latter model we make additional assumption about barotropicity of the 
process to obtain a system having an exact solution in the form of the first 
terms of the Taylor expansion in the center of vortex. Thus, we get a system of 
two equations. We could obtain the same result by assuming that the process 
is isochoric, i.e. the pressure is proportional to the temperature. We study the 
influence of the linear friction on the behavior of vortex. 

We show that in the case of atmosphere, where the vertical motion satisfies 
the hydrostatic balance and the horizontal and vertical motion are somewhat 
separated, the vortex behavior can be described by the same system of ODEs 
as the vortex in the 2D barotropic model. The vertical coordinate can be 
considered as a parameter. 

First of all, we consider a steady axisymmeric vortex in the friction-free case 
and show that in contrast to the incompressible case, there exist some param¬ 
eters such that the vortex is unstable with respect to small perturbations of 
initial symmetry in the class of solution with a linear profile of velocity. 

Further, we prove that the vortex always loses the stability when the surface 
friction coefficient is constant. We note that the flutter instability (a blowing- 
up vibrational motion) can be induced by dry friction in mechanical systems 
which would be stable without frictional forces (e.g. [26]). 

Finally, we study an axisymmetric vortex in the case of constant surface fric¬ 
tion both analytically and numerically and show that the complicated features 
of the vortex (e.g. a sudden decay and further intensification) can be explained 
by our model. Moreover, we have shown that the phenomenon of interaction 
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of the tropical typhoon with land can be quite realistically explained by a 
relatively simple ODE system, where the effect of topography is modeled by 
variable surface friction coefficient. 

The aim of our study is to show that many interesting effect of the atmospher¬ 
ical vortex motion can be qualitatively explained already by a very simple 
mechanical model. Of course, we do not pretend to claim that the trajectories 
of real tropical cyclones can be entirely explained by this model. There are at 
least three reasons for this. First of all, the real atmospheric vortex is localized 
and has structure (3.2) only in a vicinity of its center. The exact solution with 
linear profile of velocity gives the exact separation of local and leading fields 
only if the leading field is also linear with respect to the space variables. In this 
case the discrepancy q = 0 (see (3.3)) and the position of the center of vortex 
can be computed from the nonlinear system of ODEs (2.9),(4.1),(4.2), (3.1) 
exactly. If the velocity and pressure have a more realistic localized structure, 
then the discrepancy g ^ 0. Nevertheless, as follows from (3.4), q remains 
small until the divergence of the velocity of the local held is small. As we 
have shown in [21], for the case /i = 0, the difference between position of the 
center obtained from the ODE system and the position of the localized vortex 
obtained from direct numerical computations can be very small within several 
days. The second reason is that for big velocities the drag friction coefficient p 
depends on velocity itself (thus, the surface stress parametrization should be 
quadratic). Nevertheless, as one can easily see from (1.1), the quadratic drag 
friction together with the geostrophic condition plLU h + V hP = 0 (V# stands 
for the horizontal gradient) lead to a fast singularity formation for any realis¬ 
tic profile of velocity. Therefore, this kind of friction worsens the smoothness 
of the solution significantly. The third reason also relates to our assumption 
about the smoothness of solution. Indeed, as we have shown in Remark 2.5, 
the fact of the envelope of oscillations growing with time contradicts to the 
conservation of energy for the moving volume. This contradiction can be re¬ 
moved if we assume a formation of shock wave within the volume. Therefore, 
we can say that the ” theoretical” trajectory of vortex is close to a real one only 
for some period of time depending on initial data. The computations made for 
realistic parameters suggest that this interval is less than one day. Therefore 
the rising oscillations that one can see in the solution of the nonlinear system 
of ODE cannot develop. 

All these questions about correspondence of theoretical and real vortex can 
be solved only numerically. The influence of different kinds of friction on the 
trajectory of a localized vortex and comparisons with the trajectory computed 
from (2.9),(4.1),(4.2), (3.1) are the issues of our future researches. 
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